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H Abstract 

Heavy quarks are a unique probe to study the medium produced in ultra-relativistic heavy ion collisions. The dominant process 
of charm quark production at RHIC is believed to be initial gluon fusion which can be calculated in the perturbative QCD. The 
upper limit of FONLL calculation seems to be in good agreement with charm cross section measurements at mid-rapidity in p + p 
collisions at ^s NN = 200 GeV provided by STAR. The same measurement in Au+ Au collisions at equal energy reveals the number- 
of-binary-collisions scaling of charm cross section indicating that charm production is dominated by initial hard scatterings. In 
this article, we report the measurements of D°, D* in p + p at 0.6 GeV/c < p T < 6 GeV/c and D° in Au+Au collisions at 0.2 
GeV/c < p T < 5 GeV/c via hadronic decays D° — > K~7T + , D* + — > D°7T + — > K~7: + 7: + at mid-rapidity \y\ < 1. 
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1. Introduction 

The charm quark production is dominated by initial gluon fusion at initial hard partonic collisions and well de- 
scribed by perturbative QCD (pQCD) because of large charm quark mass (~ 1.5 GeV/c 2 ) Q. The mass is almost 
exclusively generated through its coupling to the Higgs field in the electroweak sector, while masses of (u, d, s) 
quarks are dominated by spontaneous breaking of chiral symmetry (CS) in QCD |2 ). This means that charm quarks 
remain heavy even if CS is restored, as it likely is in a QGP. Since charm is produced mainly in initial hard scatterings, 
one expects that charm production total cross section cr^ d N should scale as a function of number-of-binary-collisions 
A^bin- In addition, if charm quarks flow elliptically, there must have been enough interactions to easily thermalize light 
^ quarks. Hence, charm quarks are an ideal probe to study early dynamics in high-energy nuclear collisions. 

• i-H 

2. Measurement 



Charm cross section at mid-rapidity is calculated from open charm hadrons yields. These yields are obtained from 
the invariant mass reconstruction of open charm mesons through hadronic decays: D°(D°) -> (BR = 3.89%) 

and D*± -> D°(D )n ± (BR = 67.7%) -> ^ _ 7r + ^ ± . The identification of daughter particles was done in the STAR 
experiment (see Fig. [TJ at mid-rapidity \y\ < 1 at ^/snn = 200 GeV. 

The analysis presented herein was done using two datasets; the first one collected in 2009 with minimum bias 
trigger defined as a coincidence in Vertex Position Detectors East and West (105 million events total) and the second 
one collected in 2010 with minimum bias trigger defined as a signal in Zero Degree Calorimeter (280 million events 
total). 

Sub-detectors used in this analysis were: 
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Figure 1. The STAR detector. TPC (Time Projection Chamber) is main detector for tracking and PID (providing dE/dx, p), TOF (Time Of Flight) 
for PID improvement, and pile-up tracks removal and BEMC (Barrel Electromagnetic Calorimeter) in this analysis for pile-up tracks removal if 
TOF is not available. 



• Time Projection Chamber (TPC) providing 3D image of a particle track and ionization energy loss (dE/dx) and 
covering a large acceptance with full 2n azimuthal angle at \r]\ < 1. With the help of uniform magnetic field one 
can calculate momentum vector of a track as a function of track helix radius and magnetic induction. 

• Time-Of-Flight (TOF) detector providing time of the particle's flight from primary vertex to a channel pad and 
covering 72% in 2009 and 100% in 2010 of the whole barrel. When combined with momentum measurements 
from the TPC, this allows to separate pions from kaons up to 1.6 GeV/c l4l . 

STAR doesn't have any subsystem being able to reconstruct secondary vertex of D° or D* decay; one must calcu- 
late invariant mass of all Kn pair coming from vicinity of primary vertex. This results in a large combinatorial back- 
ground which was reconstructed via mixed-event method (Au+Au dataset), same-charge-sign, and kaon momentum- 
rotation (p+p dataset) and subtracted from invariant mass spectra of daughter particle pairs [5 ]. To reconstruct Z)*, 
one may exploit the softness of D* — > D°7t decay; combine low momentum pions with D° candidates, i.e. pairs with 
1.82 < M(Kn) < 1.9 GeV/c, and plot difference M(Knn) - M(Kn). Such value occupies the beginning of the phase 
space; hence it does not suffer from large combinatorial background making D* direct observation possible. The 
combinatorial background was reconstructed by side-band (picking D° candidates outside the D° mass region) and 
wrong-sign (picking soft pion with opposite charge) methods. The difference between these methods is the dominant 
source of systematic uncertainties for both D° and D* analyses. 

Open charm hadrons yields Y D o D * were calculated as Gaussian function areas from the invariant mass peaks fits 
shown in Figure [2] Raw counts were corrected with the reconstruction efficiency in used sub-detectors to obtain the 
correct Y D o D * . 



3. RESULTS 

Yields were calculated in 6 pr bins (2 for D°, 4 for D*) in p+p dataset and 5 pr bins (all for D°) in Au+Au dataset. 
Then the invariant charm cross section dcr^ c +v / '(2np T dp T dy) was calculated by formula 

d(T cc Y(T nsd 6t 



Inprdprdy 27rp T Ap T AyBR Nff rSLg 
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Figure 2. Left panel: D° signal in p+p 200 GeV collisions after same-sign background subtraction. Middle panel: D° signal in Au+Au 200 GeV 
collisions after mixed-event background subtraction Right panel: D* signal in p+p 200 GeV collisions. Combinatorial background is reproduced 
by the distributions from the wrong-sign (black dotted) and side-band (blue solid) methods. 



in each p T bin (<r NSD is non-single diffractive p+p inelastic cross section, /f rag< is the ratio of charm quarks hadronized 
to open charm mesons and 6j is the trigger bias correction). The charm cross section at mid rapidity dcr cc /dy was 
obtained from power-law function fit to dcr cc /(Inprdprdy) points (see Fig. [3J) as 170 ± 45(stat.)+^(sys.) pb. In 
Au+Au collisions we calculated invariant yield d 2 N/(N ey pTdpTdy). 




Figure 3. Left Panel: cc pair production cross section (symbols) as a function of pj in 200 GeV p+p collisions. Right Panel: The charm production 
cross section per nucleon-nucleon collision at mid-rapidity as a function of A^i n - 

The dcr cc /dy at mid-rapidity in Au+Au collisions was extracted, from a power-law fit as 1 86 + 22(stat.) + 30(sys.) + 
18(norm.) pb assuming that the /f rag . does not change from p+p to Au+Au collisions. The charm cross section for three 
centrality bins, 0-20%, 20-50% and 50-80%, is obtained according to the integrated yields. The charm production 
cross section per nucleon-nucleon collision at mid-rapidity as a function of A^in is shown in the right panel of Fig. 
[3] Within errors, the results are in agreement and follow the number of binary collisions scaling, which indicates 
that charm quark is produced via initial hard scatterings at early stage of the collisions at RHIC. The FONLL (orange 
band) and NLO 1 6 ] (light-blue band) uncertainties are also shown here for comparison. 

The R^ A depicted in Figure 4] was obtained via dividing D° yields in Au+Au 0-80% minbias collisions by the 
power-law fit to p+p yields scaled by A^i n , The uncertainty of the p+p power-law shape is taken into account as 

3 



/Nuclear Physics A 00 (2012) 7-g| 



4 




I 1 I 1 I 1 I 1 I 1 I 1 I 1 
• STAR Prel. 0-80%; 

o Alice 0-20% 



j i i i i i i i i i i i i i 

'0 12345678 



p r (GeV/c) 



Figure 4. D° nuclear modification factor Raa as a function of p. ALICE data taken from (7) 



systematic error. No suppression is observed at pr < 3 GeV/c. 



4. CONCLUSIONS 

Open charm hadrons (D°, D* + ) are measured in p+p and Au+Au collisions at y%N = 200 GeV at STAR. Charm 
cross sections per nucleon-nucleon collision at mid-rapidity follow the number of binary collisions scaling. The charm 
pair production cross sections per nucleon-nucleon collision at mid rapidity are measured to be dcr cc / (27rp T dpTdy) = 
170 ± 45(stat.)+^[(sys.) pb in p+p and 186 ± 22(stat.) ± 30(sys.) ± 18(norm.) pb in Au+Au minimum bias collisions. 
In the near future the STAR Heavy Flavor Tracker [ 8 ] will provide the necessary resolution to reconstruct secondary 
vertices of charm mesons, which will increase the precision of charm measurements. 
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